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INTRODUCTION 


<s> 


The  program  described  herein  was  devised  primarily  to  provide 
acquisition  data  to  a  ground  based  instrumentation  net  consisting  mainly  of 

9 

various  types  of  cameras.  The  cameras,  in  this  case,  are  employed  to 
photograph  a  flashing  strobe  light  that  is  mounted  on  each  pole  of  the 
ANNA  geodetic  satellite  that  is  magnetically  stabilized  along  its  spin  axis. 

Provisions  have  been  made  to  compute  a  satellite  ephemeris  by  a  dif¬ 
ferential  correction  procedure  as  the  satellite  is  stepped  around  the  orbit  in 
some  desired  increment  of  time. 

It  is  assumed  that  the  density  distribution  of  the  Earth  is  axially 
symmetric  and  that  the  force  field  is  represented  by  the  principle  term 
and  the  zonal  terms  2  through  4.  Provision  is  also  made  for  accommo¬ 
dating  the  parameters  of  a  model  atmosphere.  At  each  step  computations  are 
made  to  determine: 

1)  If  a  given  observing  site  is  in  darkness  (elevation  angle  to  the  Sun 
less  than  some  desired  €). 

2)  If  the  elevation  angle  from  the  observing  site  to  the  satellite  is 
positive  or  greater  than  some  desired  e'. 

3)  The  components  of  the  magnetic  field  (North,  East,  vertical, 
horizontal  and  total  field)  and  the  dip  and  declination. 

4)  Which  strobe,  if  either,  is  visible  to  an  observing  site. 

5)  If  the  recorded  image  size  on  the  photographic  plate  will  be  larger 
than  some  e". 

If  the  above  conditions  are  satisfied,  the  program  continues  and 
computes  the  following  additional  information  for  each  observing  site: 

1)  Time  of  the  observation. 

2)  The  azimuth  and  range  and  the  topocentric  hour  angle  and  declina¬ 
tion  to  the  satellite. 

3)  The  latitude  and  longitude  of  the  sub -satellite  point. 
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4) 


The  angle  between  the  observer-satellite  vector  and  the  observer- 
moon  vector. 

5)  If  the  satellite  is  illuminated  by  the  Sun. 

6)  The  angle  between  the  observer-satellite  vector  and  the  center  of 
the  light  cone. 

Other  features  of  the  program  are: 

1)  The  selection  of  the  most  valuable  observations  for  a  station  to 
make  from  series  of  possible  observations.  This  is  achieved  by 
considering  the  azimuths  at  which  the  station  has  previously 
recorded  data.  On  this  basis,  a  final  selection  of  flash  times  and 
the  associated  "look  angles"  from  each  site  is  made. 

2)  The  preparation  of  teletype  messages  (see  Appendix  1). 

3)  The  capability  to  compute  acquisition  data  for  observers  not 
concerned  with  the  light  (range  or  range -rate  stations). 

4)  A  station  designated  as  a  "  share  station"  will  not  have 
flashes  scheduled  specifically  for  it.  However,  acquisition  data 
v/ill  be  computed  for  any  scheduled  flashes  that  may  be 
observed  by  the  "  share  station.  " 
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II 


EPHEMERIS 


ft 

1 


The  standard  elements  (Figure  1) 

a  =  semi  major  axis 
e  =  eccentricity 
i  =  inclination 

S3  =  right  ascension  of  the  ascending  node 
<>>  =  argument  of  perigee 
M  =  mean  anomaly 
v  =  true  anomaly 

u  =  «  +  v  ,  the  argument  of  latitude 

are  used  in  the  computation  of  the  ephemeris.  Time  is  considered  the  independ 
ent  variable  and ,  as  the  satellite  is  stepped  around  the  orbit,  computations  are 
made  to  determine  the  perturbative  effects  of  the  Earth'  s  oblateness  and  drag 
on  the  satellite' s  position. 
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Figure 


The  radial,  transverse,  and  normal  perturbative  accelerations  (see 
2)  due  to  oblateness  are: 

•  • 


35 


+  —  cos  4u  sin 
48 


Hi 


3  =  “ij  -  sin2i  sin  2u  +  sini  [(1  -  |  sin2i)  cos  u  sin2icos  3u] 
-  sin2i£sin2u  (1  -  ■jr  sin2i)  +  —  sin2i  sin  4u  J  | 


•  \ 


rb 


B  r2 


2J  .  . 

—  sini  cos 
r  2 


i  sinu  +  ~s  cosi^l  -  sin2  i  (1  -  cos2u)J 


K 


-  —4  sini  cos  l 


i  £2  sinu  (1  -  ^  sin2i)  +  ^  sin2i  sin  3uJ  | 
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The  radial,  transverse  and  normal  perturbative  accelerations  due  to 
drag  are: 

r_  =  —  p  a  e  V  sin  £  £ 

D  m 

where: 


E  •- 


1  -  e  cos  E 


p  =  air  density 


A 

rv_ 

s  —  — — 

D 

m 

•  \ 

A 

rb 

D 

’  m  P 

where: 

pa 


(1  -  e2)2  V  [l  -  d(l  -  e  cos  E)2/(l  -  e2)]  E 


v  -(S)‘(K 

■(?) 


e  cos  E 


e  cos  E 


)*  [‘  - 


I 

(1  -  e2)2  cos  i 


cj  =  rotational  rate  of  the  Earth 
s 

Expressions  for  the  perturbative  effect  on  the  elements. 

2  n 


\ 

a 


3  a  n 


where 


n'  =  3  L£ 

n  "  1  '  e’Ui 


p  =  a(l  -  e2) 


ep  sin  v  + 


£i_  £] 

VJTp”  r  J 


e  =  sin  v\  +  — —  \(—  +  l\  cos  v  +  el 

^p-\r  /  ^]\T  /  J 


dJ 

Viip" 


,\  r  b 
x  =  -  cos  u 
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\ 

n 


r*b 

n/m^p 


sinu 

sini 


\  '  ' 
w  =  u  -  V 

\  V 

u  =  -  £1  cos  i 


\ 

E 


sin  E  - 


•  S  A 

r  rv 

'JiTa 


Perigee  height,  h^  =  a(l  -  e)  -  Re^ 

Radial  distance,  r  =  a{l  -  e  cos  E) 

The  difference  between  the  height  of  the  satellite  above  the  Earth  (hM)  at  any 
time  and  the  height  of  the  satellite  above  the  Earth  at  perigee  /  h  \ 


h"  -  h  =  (Re  -  Re,„  )  +  ae[l  -  cose] 

ir  ir  h  L  J 


The  air  density,  p  =  p  e 

TT 


-Kdh 


where 


K  =  i  = 


1 


h  scale  height 


p  =  air  density  at  perigee 

TT 

The  total  perturbative  effect  of  bulge  and  drag  on  the  elements  at  any  time 
a  =  a’  +  (aD  +  aB)  dt 
e  =  e1  +  (eD  +  eB)  dt 
i  =  i1  +  (iD  +  ig)  dt 
=  Q'  +  («D  +  ^  dt 
w  =  <*>'  +  («D  +  wg)  dt 
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E  =  E'  +  (E  +  ED  +  EB)dt 

u  =  u'  +  (u  +  u_  +  u_)dt 
JJ  D 

where: 

E  =  - 2 _ 

1  -  e  cos  E 

•s/1  -  e2 

u  =  1 - =:  E 

1  -  e  co  s  E 

Dots  denote  the  two  body  changes. 

Primes  denote  the  value  at  the  previous  step. 
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Ill 


ACQUISITION  DATA 


At  each  step  values  of  a,  e,  i,  fi  f  u  and  v  are  computed  after  which  the 
program  enters  the  11  acquisition  data"  phase.  The  computations  for  this  phase 
break  into  two  groups:  those  which  are  made  only  once  and  those  which  are 
made  once  for  each  observation  station. 


When  a  coordinate  system  is  not  specified  in  what  follows,  an  "Earth- 
fixed"  system  is  to  be  assumed,  i.  e.  ,.  i_and  j_in  the  equatorial  plane  with  i 
pointed  toward  Greenwich,  k  =  i  x  j_. 


Group  1:  Computations  made  only  once 

1)  Compute  X  =  west  longitude  of  node. 


a)  =  right  ascension  of  Greenwich. 

=  e  _  +  (t  - 1  y  e,  +  (t  - 1  )  e* 

Go  °  days  °fract 


b)  X  =  0_  -  d 
n  G 


-  Tf  <  X.  S.  17 
n 


2)  Compute  X  ,  west  longitude  of  sub-satellite  poi^jt. 


X  =  X  -  AX 
s  n 


-  it<  X  <  IT 

s 


3)  Compute  4>s>  the  geodetic  latitude  of  sub- satellite  point. 
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6  =  declination  of  satellite 


sini  sinu 

tan  6  =  "  ■■  -  —  ■■ 

VI  -  sin4 i  sin2  u 


tan  <$> 

s 


— ^7  tan  6;  e  2,  spheroid  eccentricity  =  .  0067686579 

1  -  e2  * 


4)  Compute  r,  distance  from  center  of  Earth  to  satellite, 
r  =  a(l  -  e  cos  E) 


Group  2:  Computations  made  for  each  station  at  each  time  step 


5)  Elevation  angle  to  the  Sun 

Elev  =  sin-1  [  sin  <f>  sinS0  +  cos  <j>  COS6(9  cos  p] 
where  4>  =  observers  latitude 

6_  =  declination  of  the  Sun 

© 

P  =  LST  -  a0 

LST  =  local  sidereal  time 
a0  =  right  ascension  of  the  Sun 


where 

a)  a©  =  -  C16  sin  2fo ;  60  =  C17  sin  a@ 

b)  L0  =  longitude  of  Sun  at  tQ  (Jan  0.  0  of  year  of  interest) 


c)  =  L0  +  C2  (t  -  40)  +  c»  sin  <c2  [ 1  -  *,1  *  C14> 


d) 


a 


6)  If  the  satellite  is  illuminated  by  the  Sun,  the  angle  I  must  be  positive  (Figure  4a). 
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p 


A  =  sin'1  (^) 

d  =  cos'1  (sin  (j>e  sincj>0+  cos  <|>»cos  4>#cos  p  ) 
where:  0  denotes  satellite 


§  =  |rab  -  RAe| 

7)  Magnetic  dip  and  declination  are  computed  by  evaluating  the 

gradient  of  the  magnetic  potential  at  altitude  (details  in  Appendix  II) 
to  obtain  the  x,  y  and  z  components  of  the  magnetic  field.  These 
quantities  permit  the  computation  of  (Figure  5). 


a)  Horizontal  component  of  Earth1  s  magnetic  field  (H) 

b)  Total  field  vector  (F) 

c)  Vertical  component  of  Earth’s  magnetic  field  (V) 

d)  Horizontal  component  of  Earth’s  magnetic  field  (H) 

e)  Total  field  vector  (F) 


Magnetic  dip,  am 


Magnetic  declination,  6  m 


cos"1 


plus  Easterly 


8) 


The  light  angle,  9,  (Figure  6)  is  the  angle  between  the  station- 
satellite  vector  and  a  unit  vector  in  the  direction  of  the  center 
of  the  light  cone.  Orientation  of  a  unit  vector  in  direction  of 
the  North  axis  of  the  light  cone  as  a  function  if  its  latitude  (<j>), 
longitude  (X),  magnitude  dip  (am)  and  magnetic  declination  (6m)>  is 


^  =  iicosa  sin6  sinX  -  cos  X  [cos6  cosa  sin<b  +  sina  cosd)]? 
— )  mm  mm  m  ) 

+  j  <cosa  sin6  cos\  +  sinXfcos  6  cosa  sin d>  +  sina  cos«b 

mm  m 


+  k  <cos  6  cosa  -  sina  sin  4>/ 
—  1  m  m  m  } 


e  = 


cos 


A 
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Define  a  unit  vector  Z  normal  to  the  spheroid  at  the  observer  (Figure  7) 


Z  =  cos  4>  cos*  jj,i_+  cos  4>  sinX^,j_+  sin  <f>  k 
A  vector  N,  in  the  meridian  and  pointing  North 
N  =  -sin  <j>  cos^  i  +  sin  4>  sin*-  i  +  cos  <|>  k 

—  E*"  -hj-*-  ““ 

E  =  Nx  Z 

u  is  a  unit  vector  from  the  station  to  the  satellite 


u  = 


TS 

ITS} 


(see  Figure  6) 


Elev  Angle  =  sin"1  (Z  •  u) 

10)  Azimuth  from  the  station  to  the  satellite 

-l  I  -x- 

Azimuth  =  cos  \  — - - 

ll)  Topocentric  hour  angle  and  declination. 

L,  a  unit  vector  from  the  station  to  the  satellite  in  inertial  coor¬ 
dinates  is 

L  =  x  i_+  y  j_+  zk 

Declination  =  sin'Vz) 

Right  ascension  =  tan 


Hour  angle  =  LST  -  RA 


Some  camera  stations  require  the  hour  angle  and  declination  of  each 
flash.  Since  the  ephemeris  is  stepped  in  time  by  an  amount  not  equal  to  the 
time  interval  between  flashes,  but  by  a  multiple  of  the  interval  between  flash 
sequences,  the  time  derivatives  of  x,  y,  and  z  are  computed  and  multiplied 
by  the  time  between  flash  intervals  to  obtain  the  new  coordinates  of  each  flash. 
With  these  new  coordinates,  the  right  ascension  and  declination  are  then 
:omputed  as  noted  above. 
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The  time  derivatives  are 


•  •  2  2 

x  s  r[l  -  sin  i  sin  u]2  cosX.  + 


\  ^  i 

r  [(-  sini  cos  i  sin2ui  +  sin2i  sinu  cos  u  u)/(l  -  sin2i  sin2u)2  cosX  ]  + 


1  \ 
r  [  (1  -  sin2i  sin2u)2  sinX/(cos2ucos2i  sin2u)  (cos  i  u  - 

\  V  ^ 

sini  sinu  cosui)  +  S2  -  u>  J 

s 

1 

y  =  -  r[  1  -  sin2i  sin2u]2  sinX 


12)  Image  size  of  the  flash 

Image  diameter  (in  microns)  =  aj  q  +  a2  q  2  +  a3  q3 

where  at,  a2,  a3  are  functions  of  the  plate  emulsion.  For  103F  emulsion 
developed  8  minutes  aj  =  7.468,  a2  =  -0.  112237,  a3  =  .0008352 

q  =  (T  B  )*  e"°‘  464  m  /8inh 

D  =  aperture  of  camera 

.«  =  distance  from  satellite  to  observer  =  |  TS  | 

T  =  transmission  of  lens 
B  =  9000  if  0  <  45° 

150(105  -  0.)  if  0  >  45° 

Am  =  atmospheric  extinction  at  the  zenith  for  the  station;  when 

an  actual  measurable  value  is  unavailable.  Am  =  1.  25  (the  • 

factor  for  moderate  haze)* 
k  *  elevation  angle 
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13)  Right  ascension,  a  ,  and  declination,  5  ,  of  the  moon  are  obtained 

m  m 

from  a  lunar  ephemeris  listed  in  memory,  a  and  6  are  given 

m  m 

for  0  hours  U.  T.  of  each  day  of  the  year.  A  six  point  interpolation 
yields  the  desired  quantities  at  any  time,  t.  For  purposes  of  this 
program  the  vector  from  the  center  of  the  Earth  to  the  Moon  is 
taken  as  being  the  same  as  the  vector  from  the  observer  to  the 
Moon.  The  angle  between  the  station-satellite  vector  and  the  station 
moon  vector,  TM  ,  is: 


Moon  Angle 


Moon  Phase,  P  ,  is  determined  from 
m 


P  =  R.  A.  -  R.  A.. 
m  m  ® 


if  P  =  0°  i  45°,  Phase  new 


m 


P  =  180°  ±  45°,  Phase  full 
m 


Otherwise,  Phase  quarter 


14)  Observation  Selection.  From  a  series  of  possible  observations, 
one  observation  (the  least  valuable  one)  is  discarded.  If 
a  station  has  made  n  good  observations  at  azimuths  A^,  A^,  A^  -  -A 
and  it  is  now  possible  to  make  an  observation  at  azimuth  B*  a 
"value"  is  assigned  to  azimuth  B  in  the  following  fashion: 


form  A  =  —  Y  '  A. 

n  1 

i  =  1 


and 


n 


P  = 


n 


(A.  -  X)Z 


i-  1 
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take  as  the  "  value"  of  the  observation 


V(B)  = 


Note  that  0  i  V  (B)  i  1,  that  V  (B)  is  zero  if  A  »  B,  and  that  it 
also  is  small  if  p  is  large  (that  is,  there  is  a  large  "  scatter"  in 
the  azimuths  already  observed).  The  values  n,  p,  and  A  are 
"updated"  as  additional  good  observations  are  made. 
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INPUT  AND  OUTPUT  FORMATS 


A.  For  Determining  Acquisition  Data  and  Selecting  Observations 

The  OGP  program  uses  logical  tapes  8  and  11  for  input.  Logical  tapes 
9,  10,  12  and  14  are  used  as  intermediate  binary  tapes. 

The  acquisition  data  and  selected  observations  are  available  as  printed 
output  on  logical  tape  5  with  data  select  zero.  With  data  select  4  the  teletype 
messages  to  APL,  NASA  and  the  1381st  Geodetic  Survey  Squadron  are  punched 
from  tape  5.  (Teletype  messages  are  also  available  as  printed  output  on 
logical  tape  6  with  data  select  zero).  Control  car^ds  necessary  for  this  comput¬ 
ation  are  input  cards  1  through  6  (Table  I  )  and  four  station  cards  for  each 
observing  site  involved  (Table  II  ).  Note  that  in  the  preparation  of  station 
cards  (Table  II)  a  station  may  be  designated  as  a  "  Share  Station.  " 

B .  ForSphemeris  Computation  Only 

If  so  desired,  the  ephemeris  portion  only  is  available  in  the  form  of 
latitude,  longitude  and  time.  In  this  instance  the  first  six  control  cards 
must  be  followed  by  four  blank  cards. 

• 
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TABLE  I 

INPUT  FORMAT  -  ELEMENT  CARDS 


Card  No. 

Card  Column 

Format 

Contents 

1 

1-6 

12A6 

Year,  month,  day  of  this  run 

7-72 

Job  heading  # 

0 

2 

1-12 

FI 2.  9 

EPOCHT  (days  and  decimals  of  day) 

13-16 

14 

IYEAR  (year  of  epoch) 

• 

17-22 

F6.  3 

CONANG  (\/Z  light  angle  in  deg.  ) 

23-28 

F6.  3 

SUNTST  (Test  for  sun  elev.  in  deg.  ) 

35-40 

F6.  3 

DELTAE  (AE  in  deg.  ) 

42-56 

F 1 5.  15 

QJ2  (J2)  • 

57-71 

F15.  15 

QJ3  (J3) 

3 

1-10 

F10.  6 

WO  (u>Q  in  deg.  ) 

12-21 

F10.  6 

RAO  ip  in  deg.  ) 

23-32 

F10.  5 

AXIS  (semi-major  axis  in  nautical  miles) 

34-43 

F 10 .  9 

ECCNO  (eccentricity) 

45-54 

F10.  8 

XINC  (inclination  in  degrees) 

56-61 

F6.  3 

QMASS  (mass  of  satellite  in 

• 

kilograms) 

63-72 

FIG.  4 

AREA  (cross  sectional  area  of  satellite  in 
square  centimeters) 

4 

1-5 

15 

NSTA  (the  number  of  stations  to  be 
considered) 

6-10 

15 

NREV  (number  of  revolutions  to  be 
considered  for  this  run) 

15-25 

F11.7 

DM1  (M  the  initial  mean  anomaly  in 
degr eesf 

5  • 

1-10 

110 

IRM  (the  maximum  number  of  flash 
sequences  to  be  allowed  per  revolution) 

1 1-20 

no 

ILM  (the  maximum  number  of  flash 
sequences  to  be  allowed  per  load) 

23-34 

F12.  9 

GDT  (integer  by  which  the  interval  between 

clock  pulses  is  multiplied  to  determine 
time  interval  at  which  acquisition  date  is 
computed;  the  elements  are  integrated 
at  one  half  this  interval) 
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40-49 

110 

INF  (initial  number  of  flash  sequences 
executed  to  date) 

50-59 

no 

INFL<  (initial  number  of  injections  executed 
to  date) 

60-71 

F12.  9 

TIMF  (epoch  of  an  even  flash  time) 

72 

11 

•  Blank  if  refined  projections  -  "one"  if 
long  range  (an  identifier  for  message  to 

NASA) 

1-12 

F 1 2.  9 

DDTT  (time  in  seconds  between  flashes 
in  a  sequence) 

13-24 

F12.  9 

TSTOP  (stop  time  in  days) 

25-28 

F  4.  1 

Minimum  angle  between  station- satellite 
vector  and  station  new  moon  vector  in 
^egr  ee  s 

29-32 

F4.  1 

Minimum  angle  between  station-satellite 
vector  and  station  quarter  moon  vector  in 
degr  ee  s 

33-36 

F  4.  1 

• 

Minimum  angle  between  station-satellite  vector 
'  and  station  full  moon  vector  in  degrees 

39-50 

F12.  9 

ELTIM  (time  in  days  at  which  acquisition 
data  computations  begin) 

53-57 

F5.  ,2 

Minimum  image  size  in  microns 

60 

11 

Blank  for  acquisition  data;  1  for  ephemeris 
computation  only 
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TABLE  II 
STATION  CARDS 


Card  No. 

Card  Column 

F  ormat 

Contents 

1 

1-24 

4A6 

Station  number  and  name  (number  in  col. 

1-4  with  leading  zeros) 

25-34 

F10.  5 

XLAT  (station  latitude  in  degrees) 

35-44 

F 10.  5 

XLONG  (Station  west  longitude  in  degrees) 

45-54 

F 1 0 .  2 

HEIGHT  (height  of  station  above  MSL^in 
feet) 

70 

11 

ft  J  It 

if  this  is  MOTS  station 

72 

11 

II  J  1 

if  this  is  camera  station 

• 

"  2" 

if  this  is  range  station 

11311 

if  this  is  range  rate  station 

2 

1 

11 

11  1 11 

if  this  station  can  only  share  flashes 

11-20 

F10.  7 

DIAM  (camera  aperture  in  mm) 

21-30 

F10.  7 

TRANS  (lens  transmission  factor) 

31-40 

F10.  7 

QM  (atmospheric  extinction  at  observer  zenith 

41-50 

F10.  7 

A1 

Constants  relative  to  film  type  and 

51-60 

F10.  7 

A2 

development  procedure 

61-70 

F10.  7 

A3 

3 

1-15 

F 1 5.  10 

C 

(Average  azimuth  of  all  flashes  observed 
from  station) 

16-30 

FI  5.  10 

D 

(Weight  factor  for  selecting  observations 

4 

Lowest  elevation  angle  (degrees)  attainable 
as  a  function  of  azimuth 

1-3 

(13) 

Elevation  Angle  0°-20#  Azm  (0*N) 

4-6 

(13) 

"  "  20-40 

7-9 

(13) 

"  "  40-60 

10-12 

(13) 

"  "  60-80 

13-15 

(13) 

"  "  80-100 

16-18 

(13) 

"  "  100-120 

19-21 

(13) 

"  "  120-140 

19 


22-24 

25-27 

28-30 

31-33 

34-36 

37-39 

40-42 

43-45 

46-48 

49-51 

52-54 


(13)  Elevation 

(13) 

(13) 

(13) 

(13) 

(13)  • 

(13) 

(13) 

(13) 

(13) 

(13) 


Angle  140-160 

"  160-180 

"  180-200 

"  200-220 

"  220-240 

"  240-260 

"  260-280 

"  280-300 

"  300-320 

"  320-340 

"  340-0 
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GLOSSARY  OF  SYMBOLS  USED 


A 

C 


14 


C15 


°16 

C17 

J 

H 

K 

f  a 
m 

n 

Re 

77 

Reh” 

V 

0. 


Cross  sectional  area  of  satellite  in  square  centimeters 

(longitude  of  perigee  of  sun- longitude  of  sun  at  January 
0.  0  for  year  of  interest) 

2  (eccentricity  of  earth'  s  orbit) 

tan^  (mean  obliquity  of  ecliptic/^) 

tan  (mean  obliquity  of  ecliptic) 

1623  x  10'6 

-6 x  10'6 

9  x  10'6 

longitude  of  sun  at  any  time 

mass  of  satellite  in  kilograms 

Mean  motion  of  satellite  (GM)  ^/a3^ 

radius  of  earth  at  perigee 

radius  of  earth  at  sub-satellite  position 

Satellite  velocity 

0°.  98564724 

360°.  98564724  * 


GM 


Satellite  latitude 


East  longitude 
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I 

Program  Constants  and 
Conversion  Factors 


I 


READ  INPUT  CARDS  AND  STORE  OUTPUT  TITLES 

EPOCHT 

= 

Epoch  of  Elements  in  Days 

IYEAR 

= 

Year  of  Epoch 

CONANG 

= 

1/2  Cone  Angle  of  Strobe  Light 

SUNTST 

= 

Maximum  Sun  Elevation  Angle 

DELTAE: 

X 

AE 

QJ2 

= 

J2 

QJ3 

= 

J3 

WO 

= 

“o 

RAO 

= 

AXIS 

=5 

aQ 

ECCNO 

= 

eO 

XINC 

>o 

QMASS 

= 

Mass  of  Satellite 

AREA 

X 

Area  of  Satellite 

NSTA 

= 

No.  of  Stations 

NREV 

= 

No.  of  Revolutions  to  be  Considered 

DM1 

= 

Mo 

IRM 

= 

Max.  No.  of  Flashes/Rev. 

ILM 

= 

Max.  No.  of  Flashes/Load 

INFL 

= 

Initial  No.  of  Flash  Loads 

TIMP 

= 

Epoch  of  an  Even  Clock  Pulse 

IP 

= 

Flag  to  Define  Long  Range  or  Refined  Predictions 

DDTT 

= 

Time  Between  Clock  Pulses 

TSTOP 

= 

Final  Epoch 

QLU1 

= 

Min.  Angle  Between  Satellite  &  New  Moon 

QLU2 

=■ 

Min.  Angle  Between  Satellite  &  Qtr.  Moon 

Q 

= 

Min.  Angle  Betwen  Satellite  &  Full  Moon 

ELTIM 

= 

Time  Look  Angles  Start 

QIMSZ 

Min.  Image  Size 

I 


Call  CONGET 
for  Annual  Constants 


To  r.2 
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.CONANG  =  0, 


SUNTST  =  0 


Is  HTST  =  0 


Write  Input  on  Output  Tape 
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Is  Mq  =  0 


^  Solve  Keplers 
Equation  for  E0 


NSTA  >100 


IXCES  =  0 


IXCES  =1 
NSTA  =  100 


READ  STATION  CARDS 

TITLE 

Station  Name 

XL  AT  * 

4> 

XLONG  = 

X 

HEIGHT  = 

ELEV 

IMOT 

Flag  for  MOTS 

IITYPE  = 

Station  Type 

DIAM 

TRANS  = 

Camera 

QM 

Characteristics 

A1  s 

at  This 

A2 

Station 

A3 

C 

Avg  Azm  of  Flashes  from  Station 

D 

Weight  Factor 

LOOK  = 

Azm  vs.  Elev.  Station 

1  SHARE  = 

Designates  Share  Station 

Start  A  New  Line  On  Output  Tape  and  Write  Stations 
to  be  Considered  for  This  Run  and  the  Input  Parameters 
Associated  with  Each  Station 
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Write  on  Output  Tape: 
"PROGRAM  GIVEN  MORE 
STATIONS  THAN  IT  CAN 
HANDLE -WILL  CONTINUEl 
WITH  FIRST  100" 


|  Initialize 

Parameters  Controlling 

Flow  of  Logic  and  Set  Revolution 

Counter  to  1 

Flag  1  =  0  NCT 

=  C 

Flag  2  =  0  M 

=  1 

Flag  3  =  0  IFIRST 

=  1 

TIME 

=  EPOCHT 

DT 

=  TT  IMP/86400 

DDT 

=  DDTT/86400 

INFR 

=  INFL 

JAR 

=  INF 

JAJ 

=  INF 

IELTI 

=  1 

CRAR 

=  0 

IR 

=  0 

IL 

=  0 

ITWO(J) 

=  0 

J 

=  1 

,  NSTA 

1 

- 

Compute  True  Anomaly  and  Argument 

of  Latitude 

Write  on  Output  Tape: 

Beginning  Revolution  No,  X 

IL 

TL(IL)  = 
SL(IL)  = 
VL(IL)  = 

IL  +  l 
TR(L) 

SR(L) 

VR(L) 

317  , 

t 

Continue 
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Compute  Time  of  1st  Flash 
of  All  Sequences  Scheduled 
for  This  Load  Time  and  Write 
Them  on  Tape  6  as  Message 
to  APL 


Call  Subroutine  Apend  to 
Write  End  of  Message  to  APL 


Cal  l  Subroutine  MH  to  Write 
Headerto  NASA  on  Tape  6 


Rewind  Tapes  9. 10, 12  and  Set  Up  Restart  j 
Elements  and  Parameters  j 

A 

= 

ARS 

AZI 

= 

AZRS 

RA 

RARS 

U 

= 

URS 

ECCN 

= 

ECRS 

XINC 

= 

XIRS 

W 

= 

WRS 

DMI 

= 

DMRS 

INF 

- 

INFR 

E 

= 

ERS 

IIME 

= 

TRS 

XSINE 

= 

SINERS 

CRAR 

= 

RCRAR 

FLAG2 

= 

0 

TLP 

= 

TIME 

IL 

= 

0 

H 

= 

HRS 

IR 

= 

0 

Jo 

VRS 

FLAG1 

= 

0 

M 

= 

MRS 

NCT 

= 

0 

Write  LHA  &  Dec.  of  All 
Flashes  of  Each  Sequence 
Scheduled  for  This  Load 
Time  on  Tape  6  as 
Message  to  NASA 


Call  Subroutine  AFHEDto 
Write  Header  for  Message 
to  APCS  on  Tape  6 


IFLAG4  =  3 


Call  Subroutine  SAO  to  Write 
Message  to  SAO  on  Tape  6 


Write  Total  Flash  Sequences  and 
Total  Loads  on  Tape  5 


Go  To 
10001 


Write  Az.  &  Elev.Ji  to 
Center  Flash  of  Each 
Sequence  Scheduled  for 
This  Load  Time  on  Tape  6 
as  Message  to  APCS 


Call  Subroutine  Flex  to  Edit  Tape  6 
and  Write  TTY  Messages  on  Tape  5 
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27 


Tl  =  999999999. 

I 

DO  206  L  =  T,  NIR 


I - j  TR(LL)--1  1 


=  1L  +  1 

TL(IL) 

=  TR(L) 

SL(IL) 

=  SR(L) 

VL(IL) 

=  VR(L) 

211  | 

CONTINUE  # 


SET  RESTART  ELEMENTS 

IR 

= 

0 

XIRS 

=  OXIRS 

ILP 

= 

IL 

WRS 

=  OWRS 

ARS 

= 

OARS 

URS 

=  OURS 

RARS 

= 

ORARS 

VRS 

=  OURS 

ECRS 

= 

OECRS 

MRS 

=  IMAS 

TRS 

= 

OTRS 

DMRS 

=  ODMRS 

RCRAR 

= 

ORCRAR 

ERS 

=  OERS 

SINERS 

= 

OSINERS 

HRS 

= 

OHRS 

AZRS 

= 

OAZRS 

J 
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29 


1 _ 

Compute  Total  Perturbative  Effect  Of  Bulge  and 
Drag  on  the  Elements  at  This  Time 


_ r  — 

Set  Flags  Telling  That  a  New  Time  Has  Been  Computed 
TNEW1  =  0 
TNEW3  =  0 
TNEW4  =  0 
OLDSIN  =  sin  (u0) 

SINUZ  =  sin(u-u0) 


* 
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Compute  Azimuth  to  Sat. 
&  Test  to  See  if  View  is 
Obstructed. Compute  Sat. 
Illumination  Angle 

■  l 

Call  Subroutine 
VAL  to  Assign  Value 
to  this  Observation 

—  I 

Compute  Topocentric  Hour 
Angle  &  Declination 
of  Each  Flash  in  the 
Sequence 

~~  1 

Write  Output  Tape  5: 
STATION 
TIME 

SAT.  ELEV.  ANGLE 
SAT.  AZIMUTH 
RANGE 

SAT.  LATITUDE 
SAT.  LONGITUDE 
IMAGE  SIZE 
LIGHT  ANGLE 
MOON  ANGLE 
SUN  ELEV. 

MAG.  DIP  &  DEL 
STA.  TYPE 


APPENDIX  I 


Three  teletype  messages  are  prepared  directly  from  magnetic  tape. 

Number  One  is  to  the  Applied  Physics  Laboratory  in  Silver  Springs, 
Maryland.  It  contains  the  following  information: 

1.  Day  number  since  January  0 

2.  Year  number  since  1900 

3.  Number  of  last  injection  into  satellite  memory  +  1 

4.  The  proper  light  to  flash  at  any  scheduled  time 

5.  Time  of  the  first  flash  in  each  sequence  in  milliseconds. 

The  format  of  the  message  is  as  follows: 

ANNA  IB  OPTICAL  ---  XXX  -,XX  -  XXXXXX 


XXXXXXXXXX 
d  e 

where  a  -  days 
b  -  year 

c  -  injection  number 

•The  times  given  are  in  milliseconds  from  0  hours  of  the  day  noted  in 
the  message. 

A  zero  in  d  denotes  a  North  flash  and  a  1  denotes  a  South  flash. 

Number  Two  is  to  NASA  Goddard|Space  Flight  Center,  Greenbelt, 

Maryland.  Contents  and  format  are  as  follows: 

_  f  g  b 

XXXXX  /  XXXXXX  AFCRL  X2  PREDICTIONS 

XXXX  /  XXXXX  /  XX£  /  XX  XX  XXX  X  XXXX  X  XXX 
'~T~  j  k  T  "m  '~rT^”o  '^p~”q  r 
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f 


where 


f  =  satellite  identification  number 

g  =  year,  day  and  month 

h  =  type  of  look  angle  (R  for  refined,  P  for  long  range) 
i  =  station  number 

j  =  sequence  number  this  day 

k  =  day  number  since  Jan.  0 

hour 

time  of  day  min 

seconds  to  nearest  1  / 1 0  second 
o  =  sign  of  Local  Hour  Angle  (1  indicates  negative,  0  indicates  positive) 
p  =  Local  Hour  Angle  to  nearest  1  / 1 0  degree 

q  =  sign  of  declination  (1  indicates  negative,  0  indicates  positive) 

r  =  declination  to  nearest  1/ 1 0  degree 

Number  Three  is  to  the  1381st  Missile  Survey  Squadron,  Orlando, 

Florida.  Contents  and  format  as  follows: 

STA.  NO.  DAY  HR  MIN  SEC  PASS  NO.  SEQ.  NO.  AZ(0°N)  ELEV  (day) 

XXXX  XXX  XX  XX  XX.  XXX  XXXXXX  XXXXX  XXX.  XX  XXX.  XX 
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APPENDIX  II 


Geomagnetic  Field  Calculations  For  a  5  Degree  Polynomial 


If  one  assumes  that  the  earth*  s  magnetic  field  at  the  point  r,8  ,  and  X 
(radial  distance,  co-latitude,  and  east, longitude)  arises  from  purely  internal 
sources,  then  the  components  of  the  field  (X,  Y,  Z)  are  given  by: 


where  Pn  are  the  Schmidt  orthogonal  functions  (using  the  associated  Legendre 
functions;  see  Chapman  and  Bartels,  Geomagnetism  -  Vol.  II;  pages  606-615); 
and  g™  and  h™  are  the  Schmidt  coefficients  of  the  geomagnetic  field;  Re  is 
the  earth*  s  radius  in  kilometers. 


The  total  magnetic  field,  F,  is  then  given  by: 

F  =  (X2  +  Y2  +  Z2)  !/2 

where  F  is  calculated  in  gauss.  To  convert  to  gamma,  the  relationship 
1  gauss  »  105  gamma  is  used. 


36 


Flow  Chart  for  Geomagnetic  Field  Computations 
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ition  the  program  selects  the  most  geometrically  .  In  addition  the  program  selects  the  most  geometrically 

ble"  observation  to  make  and  automatically  prepares  *  "valuable"  observation  to  make  and  automatically  prepares 

eletype  messages  that  are  sent  to  the  observing  sites.  •  three  teletype  messages  that  are  sent  to  the  observing  sites. 


